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Abstract 

We study transport properties of an inertial Brownian particle moving in viscous symmetric periodic structures and 
driven by an oscillating signal of two harmonic components. We analyze the influence of symmetric, antisymmetric 
and asymmetric signals on directed transport and reveal the shift symmetry of the stationary averaged velocity of the 
Brownian particle with respect to the relative phase of two components of the signal. The shift symmetry holds true 
in all regimes. 
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1. Introduction 26 

27 

Recent progress in the highly controlled fabrication 28 

of small structures opens new prospects for miniatur- 29 

ization of devices, machines, engines, etc. Processes in 30 

such systems can exhibit radically different properties 31 

than at the macroscopic level. For example at the micro- 32 

scopic scale, immanently there is a world of fluctuations 33 

which cannot be eliminated or even reduced. However, 34 

it can be exploited. A good example are biological mo- 35 

tors like kinesin or dynein which exploit thermal fluctu- se 

ations for their directed movement by the ratchet mech- 37 

anism [1]. At the microscopic or mesoscopic levels, ss 

ways and means of generation and control of particle 39 

transport are important issues for both theorists El and 40 

experimentalists |3|. In literature, there are many sug- 41 

gestions and examples how to generate a directed move- 42 

ment of particles GH4). Much more difficult problem IS 43 

related to a precise control of transport. In the paper, 44 

we study an archetype of transport in (spatially) peri- 45 

odic systems which is described by a Langevin equa- 4 e 

tion. In this modeling, we know what conditions have 47 

to be fulfilled in order to generate a directed motion of a 48 

Brownian particle. Moreover, properties of this system 49 
can be experimentally verified in a setup consisting of a 
resistively and capacitively shunted Josephson junction 

50 
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device EISIO. It is possible because the underlying dy- 
namics can conveniently be described by an equivalent 
equation of motion in the Stewart-McCumber model 
(U |3 [10| Ull. In our previous papers |6), we have 
studied the system driven by a time-periodic force G(t) 
which is the simplest harmonic signal G(t) = A cos(Q0 
(or G{t) = A sin(Q0)> where A and £1 are the amplitude 
and angular frequency of the signal, respectively. We 
have shown that, when additionally a constant force F 
is applied, anomalous transport in experimentally wide 
regimes can be observed: absolute negative mobility 
near zero value of F ( a linear response regime), nega- 
tive mobility in the nonlinear response regime and neg- 
ative differential mobility. In this paper we extend the 
analysis by considering the biharmonic driving. How- 
ever, we assume that the constant force F = 0. 

The paper is organized as follows. In Sec. 2, we 
present the Langevin equation determining dynamics of 
the Brownian particle in presence of ^-correlated ther- 
mal fluctuations. Next, in Sec. 3, we address the prob- 
lem of influence of the second harmonics on transport 
of the Brownian particle. In the parameter space, we re- 
veal reach transport behavior. Sec. 4 provides summary 
and some conclusions. 



2. Langevin dynamics 

We study the motion of a classical particle of mass 
m moving in the periodic, symmetric one-dimensional 
potential V(x) = A V sin(2nx/L) of the period L and a 
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barrier height 2AV. The particle is driven by an unbi- 
ased time-periodic biharmonic force 



G(t) = A[sin(Of) + 6 sin(2Qf + 0)], 



(1) 



where e is the ratio of the second harmonic amplitude 
to the fundamental amplitude A and the relative phase 
determines the time symmetry of the system. Addition- 
ally, the particle is subjected to the thermal noise. Dy- 
namics of a such defined Brownian motor is governed 
by the Langevin equation for the coordinate x = x(t) of 
the Brownian particle which has the form lTT2l 



mx + yx = -V\x) + Git) + ^2ykT f (f), 



(2) 



where the dot denotes a differentiation with respect to 
time and prime denotes a differentiation with respect to 
the argument of the potential Vix). The parameter y is 
the friction coefficient, T denotes temperature, and k is 
the Boltzmann constant. Thermal fluctuations are mod- 
eled by the zero-mean Gaussian white noise £(0 with 
the correlation function (^(t)^(s)) = Sit - s). 

We introduce dimensionless variables. The natural 
length scale is determined by the period L of the poten- 
tial Vix). The dynamics possesses several time scales. 
We define the characteristic time To determined from 
the Newton equation, mx = -V'(x), by inserting char- 
acteristic quantities, namely, mL/Tq = AV/L; hence 
mL 2 /AV. The dimensionless variables thus read: 
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The dimensionless Langevin dynamics consequently 
assumes the form 



I + yx = -V\x) + git) + ^|2fD~ |(f), (4) 

so where 

si • the re-scaled friction coefficient y = iy/m)ro is the 

82 ratio of the two characteristic time scales, To and 

ss the relaxation time scale of the velocity degree of 

84 freedom, i.e., tl = m/y, 

ss • the re- scaled potential 

Vix) = Vix)/AV = sin(27rJc) (5) 

se assumes the period 1 and the barrier height 2, 

s? • the scaled external time-periodic force 

ss git) = a[sin(W) + s sin(26L>? + 0)] (6) 

89 where the signal has the re- scaled amplitudes a = 

90 AL/AV and s = e/AV and the dimensionless an- 

91 gular frequencies co = Qtq, 



92 
93 
94 
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• the re-scaled, zero-mean Gaussian white noise 
forces obey = 6(t-s) with a re- scaled 

noise intensity Do = kT/AV. 

In the following, mostly for the sake of simplicity, we 
shall use only dimensionless variables and shall omit 
the "hat"-notation in all quantities. 

Transport properties in systems driven by this type of 
external stimulus have been theoretically studied mainly 
in the overdamped regime (131 [141 Q2I, for moderate 
damping [16], both experimentally and theoretically for 
cold atoms in the optical lattices ||T7l [T8j [T91, and for 
driven Josephson junctions 



104 3. Influence of the second harmonic of the driving 

105 From the symmetry considerations it follows that the 

106 long-time averaged velocity v of the Brownian motor is 

107 equal to zero if it is driven only by one harmonic, i.e. 

108 when s - in Eq. ([6]). In order to generate a directed 

109 motion of the motor, one has to include the second har- 

110 monic. Therefore we pose here the question: what is the 
m influence of the second component is ± 0) of the exter- 

112 nal force git) on transport properties of the Brownian 

113 particle described by Eq. ([?]). 

114 Nonlinearity and three-dimensional phase space 

115 ix,y = x,z = cot) make the system (|4]) possible to be- 
ne have chaotically in the deterministic case (Dq = 0). 
117 Many features depend strongly on the shape of basins 
us of attraction. If we however plug the temperature on, it 

119 is very likely that we destroy the present scene of attrac- 

120 tors and release the possibility for the system to proceed 

121 not only with attractors but more importantly with the 

122 deterministic unstable orbits. This situation is extremely 

123 complicated and can change from point to point in the 

124 five-dimensional parameter space {y, co, a, s,Dq}. It is 

125 almost impossible to find all features for such a system; 

126 therefore the goal of this work is focused only on the 

127 generic influence of the biharmonicity parameter s. In 

128 fact, one is able to tangle the picture even more by set- 

129 ting the frequency of the second harmonic in git) free, 

130 but authors feel that this is unnecessarily in this very 

131 work. 

132 In the following we will fix the dimensionless temper- 
as ature to the value Dq = 0.001 and focus on the stochas- 

134 tic (not deterministic) properties. 

135 3.1. Numerical experiment 

136 In order to establish the influence of the second 

137 harmonic of the driving force on transport properties 
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138 we have carried out comprehensive numerical simula- 

139 tions. We have employed Stochastic Runge-Kutta al- 
uo gorithm of the 2 nd order with the time step of [10~ 3 
hi I0~ 4 ](2n/oj). All numerical calculations have been per- 
U2 formed using CUDA environment on desktop comput- 
us ing processor NVIDIA Tesla C1060. This gave us a 
144 possibility to speed the numerical calculations up to 
us few hundreds times more than on typical modern CPUs, 
us More details on this very efficient method can be found 
14? in the work l2Tll . 

us We focus on the asymptotic current or long-time av- 

149 eraged velocity v of the Brownian particle. Averaging 

150 was performed over 10 3 - 10 6 different realizations and 

151 over one period of the external driving force T - In/o). 

152 We choose all initial positions and velocities to be uni- 

153 formly distributed over one potential period [0, 1] and 

154 the interval v e [-2, 2], respectively. 

155 3.2. Role of symmetry in time domain 

156 Properties of the time dependent driving force g(t) 

157 in Eq. ^ determine whether the Brownian particle is 

158 transported in the long-time regime, i.e. whether v = 

159 or v ^ 0. We can distinguish two special cases of the 

160 force g(t) . 

161 (i) The first case is when there is such to that g(to + t) = 

162 g(to - 1). It means that the driving is symmetric or invari- 

163 ant under the time-inversion transformation, see solid 

164 and dotted lines in Fig. [T] 

165 (ii) The second case is when there is such t\ that g(t\ + 

166 t) = -g(t\ - t). This is the case of the antisymmetric 

167 driving, see dashed and dotted-dashed lines in Fig. [T] 

168 As a consequence, in the symmetric case (i), the sta- 

169 tionary average velocity tends to zero when the friction 

170 coefficient y tends to zero: v — > when y — > 0; if 

171 then generically v 0. It is illustrated in Fig. [2] 

172 for = 7i 12, 3n/2. In the asymmetric case (ii), the sta- 

173 tionary average velocity tends to zero when the friction 

174 coefficient y tends to infinity (the overdamped regime): 

175 v — > when y — > oo; if y < oo then generically v ^ 0. 

176 Let us note that contrary to the symmetric driving, for 

177 y — > the velocity v ± 0, cf. Fig. [2] So, it means that 

178 the transport is generated by deterministic dynamics. 

179 We consider the case of the symmetric driving with 

180 = 7i /2 for the biharmonicity s = 0.5 (see Fig. [T]) and 

181 study the role of dissipation characterized by the friction 

182 coefficient y. This is the case when for y = the station- 

183 ary average velocity v = 0. When the friction coefficient 

184 increases starting out from zero, the average velocity be- 

185 comes non-zero as is illustrated in Fig. [2] The average 

186 velocity as a function of y displays non-monotonic de- 

187 pendence exhibiting maxima and minima. Moreover, it 

188 passes through zero and the current reversal phenomena 



g(t) =sin((j£) +esin(2a;t - 
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Figure 1: (color online) Dimensionless external ac driving g{t) = 
a[sin(ojt) + ssin(2ojt + 0)] for the fundamental amplitude a = 1, 
four different relative amplitudes of the second harmonics: s = 
0.1,0.5, 1,2 and selected values of the relative phase: = (blue 
dashed), n/2 (black solid), n (red dotted-dashed) and 3n/2 (black dot- 
ted). For arbitrary values of a and s, the ac driving possesses the time 
reflection symmetry for = n/2 and = 3n/2. For = and = n 
the driving is antisymmetric. For other values of the relative phase 
the driving is asymmetric. 
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189 can be detected. Because for y = the velocity v = 

190 and for y ^ generically the velocity v ^ 0, this case is 

191 called the dissipation-induced symmetry breaking [22]: 

192 the coupling to thermal bath is enough to break the time 

193 inversion symmetry. We note that for a fixed damping 

194 y, the average velocity for the phase = 3n/2 takes 

195 exactly the opposite sign to the case = n/2. 

196 Now, let us consider the antisymmetric case = 0. 

197 For y = 0, the velocity v ± 0. The weak dissipa- 

198 tion diminishes the stationary velocity in comparison 

199 to the dissipationless case. The dependence v(y) is 

200 also non-monotonic with minima and maxima. As in 

201 the symmetric case, the case with the phase = n 

202 can be obtained from the case = by the relation 

V (0 = 7T) = -V(0 = 0). 
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Figure 2: (color online) The stationary average velocity v as a func- 
tion of the friction coefficient y is depicted for selected values of the 
relative phase 0. For = n/2, 3n/2 the driving is symmetric while for 
= 0, n it is antisymmetric. Other parameters read: a = 4.2, oj = 4.9, 
e = 0.5 and D = 0.001. 

203 

204 3.3. Arbitrary shape of driving 

205 In previous subsection we focused on specific values 

206 of the phase. Here we present the numerical investiga- 

207 tion of the 3D parameter space {0, s, y). For phases dif- 

208 ferent than just mentioned above, we reveal also asym- 

209 metric external biharmonic signals. In the Fig. |3j the 

210 average velocity is presented in color plots for four dif- 

211 ferent damping constants y = 0.01,0.1,0.9,2 (panels 

212 a-d respectively) and additionally for the overdamped 

213 limit (panel e). On the abscissa we vary the amplitude s 

214 of the second component of the signal g(t) and on ordi- 

215 nate we present phase e [0, 2n]. Light colors denote 

216 positive average velocity. Color becomes darker for val- 

217 ues of v close to zero and eventually turn to dark-gray 

218 and black for negative valued average velocities. 

219 For the weak friction the average velocity has reflec- 

220 tion symmetry v(n + 0) = v(n - 0) as we would expect 



(a) 7=0.010000 




Figure 3: (color online) Influence of the second harmonic of the exter- 
nal force g(t) on transport properties of the system |4j. Dependence 
of the drift velocity on both the relative amplitude s (horizontal axis) 
and the relative phase (vertical axis) is depicted for various damping 
constants y = 0.01, 0.1, 0.9, 2.0 and for the overdamped limit (top to 
bottom). Other parameters are: a = 4.2, Do = 0.001 and oj = 4.9. 
Black dotted lines are plotted on all panels showing the phases for 
which the driving force g(t) possesses the reflection symmetry t — » -/, 
i.e. for = 7r/2 and 3n/2. 



4 



221 from a system prepared very close to the limit of the 

222 frictionless or Hamiltonian systems, because then the 

223 relation v « sin(0 + n/2) is quite well satisfied |23|. 

224 We plotted black dotted lines on each panel to guide the 

225 reader to the point where the driving force git) possesses 

226 the reflection symmetry, i.e., for = n/2 and 3n/2. 

227 As we increase the friction coefficient system loses 

228 its previous symmetry and becomes non- symmetric as 

229 one can easily see on panels (b) and (c). In other words 

230 - in the situation where both characteristic times in the 

231 system r r and the period T of the driving take more or 

232 less the same value, the battle between periodic stimu- 

233 lation and damping (not strong enough to suppress the 

234 driving influence quickly with possibility of additional 

235 energy cumulation) causes the whole irregular dynam- 

236 ics as seen on the central panel (c) of Fig. [3] If we, 

237 however, analyze situation with strong damping the pic- 

238 ture again gains the symmetry but now of a different 

239 kind, i.e. v(n + 0) = -v(n - 0), cf. panel (e) in Fig. [3] 

240 The close inspection of all case presented in Fig. [2] and 

241 Fig. [3] leads to the important conclusion that for a fixed 

242 set of all parameters, there is the shift-symmetry of the 

243 stationary velocity with respect to the phase, i.e., 

V(0) = -V(0 + 7T). (7) 

244 This relation is a particular case of a more general rela- 

245 tion 

v(s) = -v(s) (8) 

246 which follows from the symmetry considerations. One 

247 can note that the transformation — > + n is equiva- 

248 lent to the transformation s — > -s. The same relation 

249 holds true if, instead of the second harmonics, we apply 

250 a constant force F. Then of course v{-F) = -v(F) O. 

251 Remember that for any set of parameters the stationary 

252 average velocity v = when F = or s = 

253 3.4. Controlling transport by symmetric signals 

254 We analyze the case when the external driving is sym- 

255 metric. We set the phase of the second harmonics to 

256 = 7i 1 2 (see black solid curves in Fig. [T]). We check the 

257 system response to the signal against the relative am- 

258 plitude of the second harmonics s for the range start- 

259 ing from and ending at the value higher than dou- 

260 bled base driving amplitude a. In Fig. [4] these char- 

261 acteristics are plotted for selected driving frequencies 

262 oj = 0.1,3,4,4.9. From numerical analysis it follows 

263 that the average velocity changes its sign by varying the 

264 parameter s for all inspected frequencies of the exter- 

265 nal driving. It means that the shape of the external sig- 

266 nal can control values and direction of the net veloc- 

267 ity in the system. The current reversal can be multiple 
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Figure 4: (color online) The stationary averaged velocity vs. the rela- 
tive amplitude s of the second harmonics for four values of the angular 
frequency oj of the signal g(t). Other parameters are: a = 4.2, y = 0.9, 
<p = n/2 and D = 0.001. 

268 (oj = 0.1), akin to the situation described in (24). Keep- 

269 ing s constant at a certain level usually the direction of 

270 the average motion of Brownian particles changes its 

271 sign for the different values of the driving frequency 

272 (s = 0.5 or 2.0). On the contrary there are regimes 

273 within the scanned parameter space {s, oj} where regard- 

274 less the values of oj chosen the system response is qual- 

275 itatively the same. 

276 Next we explore the transport properties for the 

277 Brownian particle moving in the viscous environment 

278 with different friction coefficients. We examine the 

279 character of the system response against the signal of 

280 the different shape which we can control by tuning the 

281 parameter s (see Fig. [T]for details). There are two alter- 

282 native limits for the viscous system behavior - Hamilto- 

283 nian where system is frictionless l25l l26l l27l and over- 

284 damped where the characteristic relaxation time for the 

285 velocity r r = 1 jy is very long. Between those two pe- 

286 ripheries there is a region of moderate damping which 

287 seems to be the most intriguing (5J |6] [28l [29). It pro- 

288 vides rich spectrum of the very interesting phenomena 

289 and therefore we are going to focus on this particular 

290 domain in the following. 

291 In Fig. [5] the reflection of the impact of different 

292 shapes and strengths of external driving for the friction 

293 constant y in the range from 0.1 to 10 can be found. 

294 This means that the characteristic relaxation time passes 

295 from 10 to 0.1. If we refer this time to the second char- 

296 acteristic time of importance for system ([?]), namely the 

297 period of the external driving T = 2n/oj ^ 1.28, one can 

298 see that the point where both characteristic times are of 

299 the same order can be identified more or less in the mid- 
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Figure 5: (color online) Logarithmic dependence of the average ve- 
locity on the friction coefficient y is plotted for four relative ampli- 
tudes s of the second harmonics of the external driving g(t). Vertical 
black dashed line marks the point of the critical value of the friction 
coefficient, for which two characteristic times, relaxation time of the 
velocity r r = 1/y and period T = 2n/aj of the driving force, are equal. 
One can easily notice rich behavior of the average velocity around this 
specific value. Other parameters are: a = 4.2, oj = 4.9, = n/2 and 
D = 0.001. 



die of the chosen region of analyzed damping constants. 351 

Indeed, after examining of Fig. [5] one can easily re- 352 

veal most exciting features around essential value of the 353 

damping constant marked by the vertical dashed black 354 

line on the plot. At low friction, the average velocity is 355 

close to zero. When we, however, increase the friction 356 

coefficient to the value of around y = 0.3 the system 357 

starts to react in a different way depending on the rela- 358 

tive strength of the second harmonic s. For strengths 359 

less then or equal to 1 the current becomes positive, 360 

while for s = 2 system reacts with the opposite sign. 361 
This gives a possibility to control the transport simply 
by varying the strength of the second source of the ex- 
ternal field. When we go even further and arrive to the 
vicinity of the critical point r r = 2n/oj, the previous 
positive valued current starts to drop, crosses zero and 

becomes negative quite steeply. Surprisingly values of 362 

average velocities for all strengths higher than 0.5 pos- ses 

sess almost the same negative values just above y = 1. 364 

Additional enlargement of the friction leads to reduce ses 

of the transport possibilities of the system. It does not see 

reach zero, but decreases of several orders of magni- 367 

tude - see panel (e) on Fig. [3] for details. By setting ses 

the strength to zero we end up with the antisymmetric 369 

force and with zero current for any value of the friction 370 

constant due to the symmetry reasons. 371 



We have explored transport properties of the Brown- 
ian particles in a symmetric potential, driven by the time 
periodic biharmonic signals. We have demonstrated 
how the symmetry of driving force influences the trans- 
port features. There exists two limits: overdamped and 
frictionless. It turns out that in those two limits differ- 
ent types of symmetry exclude transport. In the friction- 
less case the system is time-reversible, thus the symmet- 
ric driving cannot distinct the direction. On the other 
hand, in the case of overdamped motion the antisym- 
metric driving leads to zero current. In all other cases, 
as the Curie principle suggests, the particle has gener- 
ally non-zero average velocity. The closer inspection 
shows that the magnitude and sign of the current has 
complex structure in the parameter space. Typically, the 
multiple current reversals occur, when one of the system 
parameters is changed. 

In this paper, thermal noise %(t) in Eq. ^ is assumed 
to be white noise of zero correlation time. In real sys- 
tems the correlation time of thermal fluctuations is never 
zero. In many situations this approximation is very well 
but there are also situations where the white-noise ap- 
proximation fails and a different treatment based e.g. on 
the generalized Langevin equation should be used 1 30 1 . 
However, it is essentially beyond the scope of the paper 
and requires separate investigations. 

Finally, let us remind that the Langevin equation 
([2]) has similar form as an equation of motion for the 
phase difference *F = *¥(t) between the macroscopic 
wave functions of the Cooper pairs on both sides of the 
Josephson junction. The quasi-classical dynamics of the 
resistively and capacitively shunted Josephson junction, 
which is well known in the literature as the Stewart- 
McCumber model [3 QTL is described by the fol- 
lowing equation 



k f — / sin¥ 



2eJ R 



2e 



= 2e Kt)+ 2e^ R 



2k B T 



(9) 



The left hand side contains three additive current contri- 
butions: a displacement current due to the capacitance 
C of the junction, a normal (Ohmic) current character- 
ized by the normal state resistance R and a Cooper pair 
tunnel current characterized by the critical current Iq. In 
the right hand side, I(i) is an external current. Ther- 
mal fluctuations of the current are taken into account 
according to the fluctuation-dissipation theorem and sat- 
isfy the Nyquist formula associated with the resistance 
R. It is an evident correspondence between two models: 
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372 the coordinate x = x ¥ - n/2, the mass m = (h/2e) 2 C, the 

373 friction coefficient y = (h/2e) 2 (l/R), the barrier height 

374 AV = (H/2e)Io and the period L - 2n. The biharmonic 

375 signal G(t) in Eq. ([I]) corresponds to the external current 

376 I(t) . The velocity v = x corresponds to the voltage V 

377 across the junction. So, all transport properties can be 

378 tested in the setup consisting of a resistively and capac- 

379 itively shunted Josephson junction device. 
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